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The thermal disordering of the d-density wave, proposed to be the origin of the pseudogap state of 
high temperature superconductors, is suggested to be the same as that of the statistical mechanical 
model known as the 6-vertex model. The low temperature phase consists of a staggered order 
parameter of circulating currents, while the disordered high temperature phase is a power-law phase 
with no order. A special feature of this transition is the complete lack of an observable specific heat 
anomaly at the transition. There is also a transition at a even higher temperature at which the 
magnitude of the order parameter collapses. These results are due to classical thermal fluctuations 
and are entirely unrelated to a quantum critical point in the ground state. The quantum mechanical 
ground state can be explored by incorporating processes that causes transitions between the vertices, 
allowing us to discuss quantum phase transition in the ground state as well as the effect of quantum 
criticality at a finite temperature as distinct from the power-law fluctuations in the classical regime. 
A generalization of the model on a triangular lattice that leads to a 20-vertex model may shed light 
on the Wigner glass picture of the metal-insulator transition in two-dimensional electron gas. The 
power-law ordered high temperature phase may be generic to a class of constrained systems and its 
relation to recent advances in the quantum dimer models is noted. 



PACS numbers: 



I. INTRODUCTION 



A motivating factor for the present paper is the 
specific suggestion that a new broken symmetry can 
explain the pseudogap phase of high temperature 
superconductorsaB. The corresponding order is a 
particle- hole condensate of "angular momentum" 2, 
termed the (i-density wave (DDW)d. There is same tell- 
tale evidence of this unusual order parameter!!! involv- 
ing circulating orbital currents arranged in a staggered 
pattern, which is directly detected as a Bragg scattering 
signal in neutron measurements. The second motivating 
factor is to place this proposal in a wider context of many 
body theory, where a strongly correlated electronic sys- 
tem can have unconventional broken symmetries in the 
ground state. In this respect, I shall briefly touch upon 
the topic of two-dimensional electron gas in Si-MOSFET 
devices!! 

The building blocks of the low energy theory corre- 
sponding to DDW are bond currents whose arrangements 
define the various order parameters that reflect particle- 
hole condensates^. The idea is clearly similar to resonat- 
ing valence bonds (RVB)Q, where the building blocks are 
valence bonds that can be described in terms of particle- 
particle condensates. These may order in the ground 
state, j-|Or they may not, in which case one has a spin 
liquidEl. In either case, strong correlation effects are be- 
lieved to play an important role. The parallel goes fur- 
ther: while attractable RVB Hamiltonian is the quantum 
dimer modelu, a tractable model for bond currents will be 
seen to be a vertex model known in statistical mechanics 
and its suitable quantum generalization. 

The actual statistical mechanics of the DDW transition 
is richer than the mean field (Hartree-Fock) picture in 
which the ordered pattern is frozen, until the magnitude 



of the circulating currents vanishes, which is the sim- 
plest possible description of the broken symmetry phase. 
While this is reasonable deep in the ordered state, it does 
not allow for fluctuations of the order parameter. 



A natural modifica.tian that I shall describe, involving 
the 6-vertex modelErlLil, leads to striking consequences 
for the pseudogap phase: (1) As the d-density wave dis- 
orders with increasing temperature, the system exits to a 
power-law phase that is not due to any underlying quan- 
tum critical points. (2) The pseudegap transition does 
not have any specific heat anomalylla. (3) Because of the 
power-law nature of the orbital current correlations of 
the disordered state, it is likely that the electron spectral 
function exhibits a cut spectrum, unlike a Fermi liquid. 
(4) As the temperature is raised, a second transition takes 
place at a higher temperature, where the local amplitude 
of the pseudogap vanishes. It is unusual, but true, that 
the two-dimensional classical statistical mechanics of the 
pseudogap state viewed as a DDW allows for a power-law 
phase with infinite correlation length in the high temper- 
ature regime, while the low temperature phase is ordered 
with a finite correlation length. 



The outline of the paper is as follows. In Sec. II, I 
introduce and discuss the 6-vertex model for the pseu- 
dogap phase, and, in Sec. Ill, I discuss its qunatum 
generalization. Section IV is a brief analysis of the re- 
lation with the quantum dimer model. In Sec. V, I in- 
troduce the 20-vertex model on a triangular lattice and 
discuss its relevance to the Wigner glass picture of the in 
metal-insulator transition. And, finally, the conclusions 
are summarized in Sec. VI. 
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II. THE 6- VERTEX MODEL FOR THE 
PSEUDOGAP PHASE 

The ordered singlet DDW state consists of a staggered 
circulating pattern of currents flowing on the square pla- 
nar CuO-lattice, as shown in Fig. |l|. The configuration 



FIG. 1: A segment of the circulating current pattern in the 
DDW ordered state. 

is a result of juxtaposing two sets of vertices, shown in 
Fig. |^, centered on the Cu atoms in a current conserving 
manner. There are only two possible choices, resulting in 



(5) (6) 

FIG. 2: Vertices 5 and 



two distinct broken symmetry states, which are related 
by time reversal and a lattice translation. 

More mathematically, the order parameter <£>q is given 

by 



(^t( k + Q )t )^ (M)) 



i — — (cos k x 



cosfc 9 )C (!) 



where a and p are the spin indices, and Q = (7r, 7t). The 
lattice constant has been set to unity, and the operator ip p 
is the electron destruction operator. The order parame- 
ter is called a density wave because it is a particle-hole 
condensate, even though what is actually modulated is 
current, not density. The order parameter <I>q is a spin 
singlet; there is also a triplet version in which it is the 
spin current that is modulated.u The order parameter 
is called a d-wave because of the internal form factor of 
the particle and the hole, which is (cosk x — cosk y ). Of 
course, on a crystalline lattice, angular momentum is not 
a good quantum number. This is the closest we can come 
to the angular momentum 2 of a d x i_ y i wave function in 
free space. Note that for particle-hole condensates there 
are no exchange requirements governing order parame- 
ter symmetries that enslave the orbital wave function to 
the spin wave function, in contrast to superconductivity 
where the condensate is of particle-particle type. Thus, 



one can have a DDW that is either a spin singlet, or a 
spin triplet. 

In a mean field (Hartree-Fock) picture, the only way 
the state can disorder is by the thermal collapse of the 
magnitude of the order parameter, which is a highly re- 
strictive mechanism because it does not allow for fluctu- 
ations. The system may disorder long before the magni- 
tude of the order parameter collapses. To build fluctu- 
ation effects, we consider the basic building blocks that 
are bond currents between the nearest neighbor sites x 
and y = x ± x, or y = x ± y, defined by 



jx, y = -ie$ Q (^ y -h. a), 



(2) 



where e = ±1 determines the direction of the current 
flow. The status of the bond currents is identical to the 
local order parameter in Ginzburg-Landau- Wilson for- 
malism in which the partition function is a sum over all 
complexions of the order parameter weighted by the ef- 
fective coarse grained action. To include quantum fluc- 
tuations, one must supplement the theory with a suitable 
dynamics, which I shall discuss in the following section. 

I shall assume that there is a regime of coupling con- 
stant and temperature such that the bond current order 
is well developed, and the fluctuations of the magnitudes 
of the bond currents can be neglected. If this is the case, 
there should be a second transition at which the mag- 
nitude of the bond current itself vanishes. This would 
imply a second pseudogap transition at a higher tem- 
perature. The existence of a specific heat anomaly at 
this upper pseudogap transition is not entirely clear, as 
it will be seen to be a transition from a power-law or- 
dered phase. However, in any case, the temperature may 
be too large to extract it from the large phonon or other 
backgrounds. I shall therefore ignore this transition al- 
together. 

It is easy to convince oneself that low lying thermal or 
quantum fluctuations can be expected to reverse a set of 
bond currents (without changing their magnitudes) pro- 
vided no sources or sinks are generated, that is, Vj = 0. 
Given that there are two incoming and two outgoing cur- 
rents at a vertex of a square lattice, there are altogether 
2J5T = 6 possible vertices. Thus, there are four additional 
allowed vertices, beyond the two shown in Fig. 0. These 
are shown below in Fig. ||. In general, each of the 6 lo- 
cal arrangements can have a distinct energy Si, but if 
we impose the restriction that there are no net external 
currents, then 



£l — £2, £3 



^4, £5 — £6- 



(3) 



The model is then unchanged by reversing all the arrows. 
For tetragonal symmetry, there are only two independent 
energy scales £1 and £5 because S\ — £3, but for generality 
we shall assume that they are distinct. 

From a Hartree-Fock analysis]!] it can be argued that e\ 
and £5 are close in energy. This is because the vertices 1 — 
4 correspond to (p x + Pj,)-density wave (PDW) states in 
the order parameter language, and the energetics of both 
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the singlet p-density wave and the DDW are controlled 
by small pair-hopping matrix elements, smaller than the 
scale of the antiferromagnetic exchange constant J. A 
singlet Pz-density wave state has the ordering 



(4) 



where Q' = (0, 7r). Note that in this case $q, is real, but 
because of the form factor sink x , a Fourier transform to 
the real space brings out the pattern of currents shown 
in vertices in Fig. |3j, when it is superposed with the cor- 
responding p a -density wave. The order parameter <I>q, is 
the closest analog on a lattice of angular momentum 1, 
p x wave function in free space. Note, as before, that it is 
the current that is modulated, not the density. 

I shall first consider thermal fluctuations and take into 
account quantum fluctuations in the following section. 
A basic assumption I shall make is that the energy of a 
given state is a simple sum of energies associated with 
the configuration at each vertex. This is a reasonable as- 
sumption because any long-ranged interaction between 
the vertices are unlikely, as the thermal smearing will 
cause any interactions mediated by the nodal quasiparti- 
cles to be exponentially decaying. Following the conven- 
tional notation, 1 shall define. 



: UJl = U>2, 



UJ 3 =UJ4, C 



W 6 , 



(5) 



-ei/T 



where the Boltzmann factors are defined by u>i 
where the Boltzmann constant ks is set to unity. The 
partition function is 



a m +n 2 ^n 3 +n 4 c n 5 +n 6 



(6) 



where is the number of vertices of type k. The sum is 
over all arrangements that fit together continuously with- 
out generating sources and sinks. The partition function 
is precisely the partition of the fa-vertex model for which 
many exact results are knownEJO. 



A. Phase diagram 

The phase diagram is shown in Fig. ^|. The regions 
I and II have orbital ferromagnetic order due to macro- 
scopic currents and correspond to a > b + c and b > a + c, 
respectively. The region III is disordered, corresponding 
to a,b,c < | (a + b + c), but the current-current cor- 
relation function exhibits a power-law. It includes the 
infinite temperature case a = b = c = 1. This entire 
high temperatuxepcegion is on the the critical line of the 
8-vertex mode£30. The region IV is the DDW phase, 
which is an orbital antiferromagnet and corresponds to 
c > a + b. The phase boundary between region III and 
IV is given by - = — - + 1, or c = a + b, which implicitly 
defines the transition temperature T c . 

Consider a given set values of a, 6, and c. As the 
temperature increases from to oo, this point follows 
a path always ending at (1,1). This path may or not 
cross a phase boundary. The path followed in the phase 
diagram when a — b is shown as an arrow. On this 
path, at T c = T*, there is a phase transition from the 
ordered DDW phase to the disordered high temperature 
phase with a power-law correlation in the current-current 
correlation function. The temperature T* is given by 



In 2 : 



(7) 



determining, phcnomcnologically, the energy difference 
As = £i - e 5 . 




FIG. 4: The phase diagram. The regions I and II are orbital 
ferromagnets (OF). The region III is the power-law phase and 
the region IV is the orbital antiferromagnet with DDW order. 
The arrow marks a path (tetragonal symmetry assumed) from 
the low temperature to the high temperature phase with the 
pseudogap transition at T c = T* . 

In region IV, the correlation length, £, i 



In 



oo , 

2*v» n u 



„4m 



(8) 



where x — e _A and — cosh A = (a 2 + b 2 — c 2 )/2ab. As 
T — ► T~, A oc (Tc-T) 1 / 2 , and it can be seen by applying 
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Poisson summation formula that 

£-! = 8e"^, A -> 0. (9) 

B. Free energy and specific heat 

The exact fpee energy of the 6- vertex model was ob- 
tained by Lieblla, and it is quite remarkable. If we denote 
the free energy density in region III by fm(T > T c ), then 
its analytic continuation to T < T c is complex, equal to 
flv(T < T c ) +«/ S mg- The singular part of the free energy 
has only an essential singularity at T = T c , which is given 
by 

/singer'- ( 10 ) 

In fact, all temperature derivatives of the free energy ex- 
ist and are identical on both sides of the transition. An 
asymptotic expansion at T c , which is_the same above and 
below the transition, is given in Ref.EiJ. 

The essential singularity implies that there are no ob- 
servable specific heat anomalies, as any derivative of the 
free energy vanishes at this infinite order transition. The 
situation is exactly the same as that of the Kosterlitz- 
Thouless phase transition of the two-dimensional XY- 
model. The transition of the 6-vertex model is inverted, 
however, as the low temperature phase has a finite cor- 
relation length, while the the high temperature phase is 
the power law phase with continuously varying critical 
exponent. The analytic part of the free energy results 
in a broad heat capacity peak above the transitiorO for 
the XY-model, but below the transition in the 6-vertex 
model! 

The specificpheat of the 6-vertex model was also dis- 
cussed in Ref.til, but it is illuminating to plot it here. 
In Fig. we show C/Nks, where N is the number of 
vertices, as a function of T/T*, where we have set a = b. 
There is no observable anomaly at T* despite the ther- 
modynamic transition at this temperature to the ordered 
DDW phase. Nonetheless, since this is an order-disorder 
transition, the high temperature entropy must disappear 
as T — » 0. At T = co, a = b = c = 1, and the enteopy 
is that of square ice, which was calculated by LiebllS to 
be S — ^f- In | . The system looses this entropy by a 
specific heat hump clustered around a definite tempera- 
ture, which is significantly below T* . The specific heat 
per vertex is also a universal function of T/T*. That it 
vanishes at high tempeartures — true for many statisti- 
cal mechanical models such as the Ising model — simply 
reflects the bounded energy spectrum on the high energy 
side and must not be taken seriously. The high energy 
spectrum is not modeled correctly by such low energy 
effective models. As a result of the boundedness of the 
energy spectrum, the entire passage from complete order 
to complete disorder takes place over a limited range of 
temperature. However, we should probably take seriously 
the fact that the DDW transition results in a specific heat 
of order ~ 0.5fcs per vertex considerably below the DDW 



C/Nk, 




0.5 1 1.5 2 

T/T* 



FIG. 5: The specific heat of the 6-vertex model as a func- 
tion of T/T* for a — b. Note that there are no observable 
anomalies at the transition at T* , and the peak of the spe- 
cific, which is entirely due to the analytic part of the free 
energy, is significantly below the transition. 



transition and may be mistaken as the pretransitional ef- 
fect as the system enters the superconducting state from 
the pseudogap phase. The specific heat peaks at the 
superconducting transition in the underdoped regime in- 
deed appear to be anomalous and very broad. The low 
temperature behavior of the specific heat is given by 

Nk B \ T ) ' V ' 



C. Floating power- law phase 

How does the coupling between the two-dimensional 
CuO-planes change the behavior of the system? Since 
the low temperature phase is ordered with a finite corre- 
lation length, small coupling between the planes will in 
general have small effect, except close to the critical re- 
gion where the correlation length diverges and a crossover 
to 3D should take place. This can be checked explicitly 
from a mean field analysis. The coupling between the 
planes due to the orbital currents is indeed very small, 
as the magnetic field at the center of a plaquette gen- 
erated by the circulating currents can be estimated to 
be of order 10 G from the magnitude of the pseudogapH. 
There is some experimental support to this facta, as the 
correlation length in the perpendicular direction does not 
extend much beyond a unit cell. 

In contrast, the effect of three-dimensional (3D) cou- 
pling in the high-temperature power-law phase is much 
more interesting. In principle, the 2D power law phase 
can be destroyed by arbitrarily weak coupling between 
the planes. However, because of the miniscule coupling 
between the planes, due to an estimated 10 G magnetic 
field of a plaquettea, the crossover to 3D order is likely to 
take place at very long length scales. Thus, this is effec- 
tively a floating phase in which the system behaves as a 
stack of decoupled 2D power-law phases. The argument, 
in principle, is the same as that of Josephson coupled 2D 
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XY-systems in a direction perpendicular to the planesS, 
although the details are different. 

Consider the current-current coupling of a given layer, 
n, with a neighboring layer (n + 1), in the continuum 
limit, which is 



f d 2 r 

H c [n] = g J —jn{r) ■ j n +i(r), 



(12) 



where s is a short distance cutoff of the order of the lattice 
spacing. It is sufficient to consider the coupling between 
the nearest neighbor planes for our purposes. We wish 
to determine the relevancy of the coupling g under the 
renormalization group transformation. Let 



(jn( r ) • j™+lWJn(0) ' jn+l(0)) 



1 

Tlx' 



(13) 



where the average is taken with respect to the fixed 
point Hamiltonians of the uncoupled planes. Since ev- 
ery point in the power-law phase corresponds to a fixed 
point with a continuously varying critical exponent, we 
merely need to determine the critical exponent x of the 
current-current correlation function. 

The exponent of the arrow-arrow correlation function 
of the 6- vertex model (or the current-current correlation 
function in the present problem) in the power-law phase 
can be obtained from the exponent of the 8- vertex model 
along the critical line, which in turn can be obtained from 
the S z -S z correlation function of the one-dimensional 
(ID) quantum XXZ Hamiltonian. The correspondence 
can also be seen by considering Trotter-Suzuki decompo- 
sition of the ID XXZ Hamiltonian, in which the arrow 
correlations along a row of the 6- vertex model (or along a 
column) correspond to correlations along the diagonal of 
the space-imaginary time latticeES. This subtlety is of no 
consequence, as we are interested in the exponent along 
a critical line at which the correlation length is infinite. 

The S z -S z correlation has a staggered part and a uni- 
form part. The exponent for the uniform part is 2, inde- 
pendent of temperature, while the exponent of the stag- 
gered part is given by 4, whereEjO 



Assuming tetragonal symmetry (a = b), 



^ 2 



1. 



(14) 



(15) 



Therefore the exponent | varies monotonically from 1 at 
T = T* to 3 at T = oo. The exponent of the uniform 
part follows from a local conservation law combined with 
the conformal invariance in two dimensional]. Clearly, 
the slowest decay will determine the 3D crossover. At 
T* , it is determined by the staggered part since its ex- 
ponent is less than the exponent of the the uniform part. 
As the temperature increases, the exponent of the stag- 
gered part increase and crosses the uniform exponent. 



We therefore define the exponent x to be either the ex- 
ponent of the uniform part, or that of the staggered part 
of the current-current correlation function, whichever 
is smaller. Under a renormalization group transforma- 
tion, for which s — > As, where A > 1 is a scale factor, 
}n{r) • jn+iW A :c j n (r) • jn+i(r) from Eq. Thus, 
for the free energy to remain invariant, it follows that 
g -> \ 2 ~ x g; note that d = 2. Therefore, at T = T* , 
where x = 1, the 3D coupling, g, is strongly relevant. As 
the temperature T — > oo, the uniform exponent [x = 2) 
takes over and g becomes marginal. 



III. QUANTUM VERTEX MODEL 

It is interesting to ask what a quantum generalization 
of this model could be. At the very least, we must allow 
the bond currents to flip as a result of quantum fluctua- 
tions. The most local move is to flip the currents around 
an elementary plaquette of the square lattice. But it is 
easy to see that such a process can create sources and 
sinks, if we do not impose any constraints. From a path 
integral description, the configuration at an intermedi- 
ate (imaginary) time slice may contain a source or a sink 
even though the states at the initial and the final time 
slices at and (3 do not. 

I start by considering the 8-vertex model written in 
terms of Ising spin variables and then incorporate quan- 
tum mechanics by introducing a transverse field on a dual 
lattice site. The advantage is that the description is in 
terms of the familiar Ising basisEU. One must then ex- 
plicitly introduce projection operators to eliminate the 
unwanted sources and sinks. 

The 8-vertex model (allowing only an even number 
of arrows into and out of each site) can be written as 
two coupled Ising models on interpenetrating square lat- 
tices with ■ an additional four spin coupling, as is well- 
knowr£3l23. The spins are situated on the dual lattice 
consisting of the sites at the center of the plaquettes of 
the original square lattice of bond currents. An arrow to 
the right (or upward) represents adjacent parallel spins; 
similarly, an arrow to the left (or downward) corresponds 
to adjacent antiparallel spins. The mapping is shown ex- 
plicitly in Figs. ^ and ^. The additional vertices shown 
in Fig. ^ are assigned energies £7 = eg, and the number 
of such vertices must satisfy nn = n§ = rid- We shall 
denote these vertices as of type d. However, ultimately, 
we must impose an additional constraint on the Hamil- 
tonian in which the vertex configurations shown in Fig. |^ 
are removed from the Hilbert space of states. 

Quantum fluctuations are incorporated by a transverse 
field, h. The resulting model, the quantum 6-vertex 
model, is then defined by the following partition func- 
tion: 



Zc 



Tr 



-PPHP 



(16) 
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FIG. 6: The mapping of the first 6 vertices to equivalent spins 
on the dual lattice. 
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FIG. 7: The mapping of the vertices involving sinks and 
sources to spins on the dual lattice. 



where the Hamiltonian H is given by 

j,k j,k 

~ J" a ^,k a j+l,k+l a j+l,k u j,k+f (17) 
j,k 

The operators a x and cr z , are the standard Pauli matri- 
ces, and P 2 = P is the projection operator that projects 
out the sources and sinks. It is easy to check that a single 
spin flip at a dual lattice site flips the bond currents on 
the surrounding plaquette, which causes transitions be- 
tween the vertices on the four corners as shown in Fig. ||. 
The previously defined vertex weights a, 6, and c are 
given by 



,{J+j'+J")/t ^ 

,(-J-J>+j")/T 
,(-J+j'-j")/T 



(18) 
(19) 
(20) 



The new vertex weight d is 

d=e ( j - j '- j "y T . 



(21) 
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FIG. 8: The flip of the central spin from H > — flips the sur- 
rounding bound currents around the plaquette. The flipped 
currents are denoted by dashed arrows. 

The classical statistical mechanics that controls the fi- 
nite temperature phase transition is unchanged because 
there are no quantum mechanical spin flips. Since the 
trace is a sum over states from which sources and sinks 
are projected out, the model is clearly the classical 6- 
vertex model — the number of vertices of type d, = 0. 
We can arrive at the same conclusion by considering an 
imaginary time path integral from a Trotter decomposi- 
tion along the imaginary time direction, t, between and 
/3, with periodic boundary condition. The constraint can 
be enforced by including a delta functional S(V ■ j(x, r)) 
at each space-time point. In the classical limit, we are left 
with a single time slice, because quantum spin flips due to 
the transverse field are punished. The partition function 
is now a sum over all spin configurations subject to the 
condition that there are no sources or sinks. Of course, 
at any finite temperature, the effective parameters of the 
model will be renormalized by quantum fluctuations at 
short distances, but that can not change the universality 
class of the finite temperature phase transition. 

At T = 0, the model is somewhat complex because of 
the projection operators. At h = 0, the ground state is 
the ordered DDW state with no quantum fluctuations, 
if the vertices c have the lowest possible energy. As the 
transverse field h increases from zero, quantum fluctua- 
tions increase and nucleates PDW vertices. It is not clear 
what the ultimate fate of the system is. There are four 
distinct possibilities: (1) a quantum phase transition to 
a gapped quantum disordered state, which I find unlikely 
because of the constraints. If this is the case, the con- 
straints have to be somehow irrelevant; (2) a transition 
to a PDW state, which I also find unlikely because there 
are no flippable plaquettes, and hence the state is disfa- 
vored from kinetic energy considerations; (3) a transition 
to a power-law phase since a local deformation may not 
be able to heal sufficiently fast. (4) There is also a re- 
mote possibility that the system does not loose its DDW 
order. In any case, the model is sufficiently complex that 
at present I cannot provide a quantitative analysis. It is 
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probably best studied by ouantum Monte Carlo methods 
involving loop algorithmstj. 

I shall assume that it is the third possibility dis- 
cussed above that is realized. If this is correct, such a 
quantum critical point is unusual from the conventional 
perspectiveEj. To illustrate this point, consider a phase 
diagram in which a system undergoes a continuous phase 
transition as a function of tuning parameters, tempera- 
ture, T, and a coupling h, such as the transverse field. 
This is shown in Fig. |9|, where the shaded region is the 
domain of classical critical fluctuations extending all the 
way to zero energy. Any influence of the quantum crit- 
ical point at h c has to be outside this domain. So, we 
may argue that, on some scale, the disordered phase at 
a finite temperature could be influenced by the quantum 
critical point. 



Disordered 




FIG. 9: Conventional phase diagram in which a finite temper- 
ature phase transition exists — the solid curve. The shaded 
region is the domain of classical critical fluctuations. The in- 
fluence of the quantum critical point at h c , at T = 0, has to 
lie outside this domain. 



The phase diagram shown in Fig. [l^, pertinent to our 
present problem, is strikingly different. Here, the exis- 
tence of a classical power-law phase of the 6-vertex model 
has entirely eliminated the quantum fluctuations, which 
are, strictly speaking, cutoff by the thermal length, while 
the classical power-law correlations have no such long- 
distance cutoff. Within the ordered phase, there is a clas- 
sical critical region, but otherwise the ordered phase is 
under the influence of the broken symmetry fixed point. 
Of course, we are assuming that the magnitude of the 
order parameter can be assumed to be fixed. If not, 
there will be a further transition at a higher tempera- 
ture where its magnitude collapses. The region above 
this transition can now be influenced by the quantum 
critical point. Nonetheless, there may be a sizable region 
of the phase diagram in which the classical power-law 
phase dominates. 



Disordered 
classical power-law phase 




FIG. 10: The phase diagram in which the entire finite temper- 
ature disordered phase is a classical power-law phase — the 
shaded region. The dashed curve is the crossover scale for the 
classical critical fluctuations within the ordered phase. The 
quantum critical regime is entirely eliminated, except pre- 
cisely at T = 0. The ordered regime is under the influence of 
the stable broken symmetry fixed point. 



IV. RELATION TO THE QUANTUM DIMER 
MODEL 



As remarked in the Introduction, there are some sim- 
ilarities with the RVB ideas, particularly with the 
cent developments involving quantum dimer modelsE 
but the distinctions are equally important. At high tem- 
peratures, the quantum dimer model on a square lattice 
maps on to— the classical dimer model with power-law 
correlationstll, which is a c = 1 (conformal charge) the- 
ory, equivalent to two massless Majorana fermions. It is 
also the fully frustrated Ising model at T = 0. 

This critical phase Wfises from the Rokhsar-Kivelson 
quantum critical pointBES. Instead, in the present case, 
the quantum vertex model maps on to the classical 6- 
vertex model, which undergoes a phase transition with 
an essential singularity from an ordered region IV to the 
power-law region III. But this power-law phase does not 
stem from a quantum critical point, instead the entire 
region is situated on the critical line of the classical 8- 
vertex model. 

Recent thrust in quantum dimer models has focused 
on the existence of a spin liquid phase, which is made 
plausible by the fact that the classical dimer model on 
a triangular lattice is non-jcritical with exponentially de- 
caying dimer correlations^. Thus, the ground state of the 
triangular case has no gapless collective excitations, only 
deconfined, gapped spinons for a finite range of param- 
eters. In contrast, we have focused on the existence of 
a power-law phase above the DDW phase and the non- 
existence of the specific heat anomaly at the transition. 
We have also argued that it is also possible to have a 
quantum phase transition in the ground state, which is 
induced by a sufficiently large transverse field mimick- 
ing the actual pseudogap transition in the groundstate. 
There are also interesting resemblances of our work with 
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the quantum and sliding ice pictures of Ref. 



V. VERTEX MODELS ON A TRIANGULAR 
LATTICE 

Vertex models may be useful in other contexts. For 
example, consider the ideas that were put forward re- 
cently regarding the metal-insulator transition in two- 
dimensional electron systems in which the transition is 
identified as a quantum phase transition in a disor- 
dered Wigner crystal stateE£lE3. The chain of reason- 
ing is as follows: extensive calculations of multiparti- 
cle exchange Ha»iltonians in pureCJCjH and disordered 
Wigner crystalsEJ suggest that there exist ground states 
with well-defined bond currents. These bond currents 
can be represented by vertex models on a triangular lat- 
tice. It is notable that recent studies have revealed that 
a metal-insulator transition could occur in such models. 
A protypical example is provided by a two dimensional 
tight-binding Hamiltonian in a random magnetic field 
and on-site disorderESES. It is in this context that it 
is useful to consider a vertex model of bond currents 
on a triangular lattice and study its quantum dynam- 
ics. Although even in- then-classical limit, such models 
are unsolved in generalt2lE3E3, some exact results, known 
for certain special cases, are sufficiently encouraging and 
similar to the 6-vertex model to pursue further. 

If we restrict ourselves to 3 incoming currents and 3 
outgoing currents at a site of a triangular lattice, this 
generates a 6!/3!3! = 20 vertex model. If the vertices 
which differ only by rotation and reflection are treated 
alike, they may be classified as follows: (i) 6 vertices in 
which the incoming arrows are adjacent, (ii) 2 vertices 
in which the incoming and the outgoing arrows alternate 
round the vertex, (iii) 12 vertices containing two incom- 
ing arrows directly opposite each other. The examples 
are shown in Fig. O. 




(D (2) 




(3) (3) 
FIG. 11: Examples of three types of allowed vertices. 



Unlike the 8- vertex, or the 6-vertex, model on a square 
lattice, the solution to the 20-vertex model is not known 
in general. It is only solved for certain values of the vertex 
configuration probabilities. As with the 6-vertex model, 
this model undergoes a phase transition to a power-law 
ordered phase, with an essential singularity in the free en- 
ergy, from an ordered state with finite correlation length 
in which the arrows alternate along a row. Since the exact 
solution is not known in general, it wauld be extremely 
useful explore this model numerically^. 

The quantum dynamics of this model can once again 
be implemented by going over to a Ising spin model on 
the dual lattice with a transverse field at a site. It would 
be interesting to explore this model, which could pro- 
vide us with a novel perspective on the old problem of a 
Wigner crystal, which, in the conventional picture, is an 
ordered state with at most small zero-point oscillatinn.s of 
an electron about a lattice site; see, however, Refs.ol23. 



VI. CONCLUSIONS 

I have argued that the natural generalization of the 
Hartree-Fock model of the DDW are certain vertex mod- 
els involving bond currents. The finite temperature tran- 
sition from the ordered DDW phase to the disordered 
phase belongs to the universality class of the 6-vertex 
model and is striking because of the absence of any ob- 
servable specific heat anomalies. I have also shown that 
the disordered phase is a power-law phase of the bond 
current correlations. These power-law correlations do 
not arise from any quantum critical point, but this en- 
tire phase is poised on the critical line of another hidden 
model, the 8-vertex model. The origin of the power-law 
correlations is the constraints that the configurations of 
the 6-vertex model must satisfy. Thus, a local deforma- 
tion cannot heal within a finite correlation length. In 
this respect, there is a strong similarity with the RVB 
dimer models on a square lattice, where the constraints 
in the dimer model are once again responsible for the 
high temperature power-law correlations. 

It is important to note that the mere knowledge of the 
nature of the order parameter does not necessarily deter- 
mine the universality class of the phase transition. For 
this, it is necessary to know also the nature of the pos- 
sible excitations that can be thermally populated. Thus, 
labelling the DDW transition as an Ising transition is an 
oversimplified view. 

I was also able to construct a quantum generaliza- 
tion involving a 8-vertex model in a transverse field, but 
with projection operators that ultimately project out the 
sources and the sinks of the 8-vertex model. The model 
is sufficiently complex to present any in-depth discussion 
at this time, but it is not difficult to motivate a quan- 
tum phase transition as a function of the transverse field. 
Nonetheless, it is remarkable that the effect of the quan- 
tum critical point is completely swamped by the classical 
power-law correlations at any finite temperature. This is 



9 



a strict departure.-£rom the conventional theory of quan- 
tum critical pointed. 

I have also speculated on the possible application of a 
similar vertex-model to the Wigner crystal and its im- 
plications to the metal-insulator transition in the two- 
dimensional electron systems. Similarly, a zoo of new 
sets of density wave theoriescl can be usefully formulated 
in terms of vertex models opening up a whole class of 
questions in the theory of quantum statistical mechan- 
ics. Some of these would .correspond to modulations of 
the bond kinetic energyO'L3. In fact, it is quite possi- 



ble that certain regions of, 
described by these models! 



le cuprate phase diagram are 
but the aim here was to con- 



sider only those models that are natural generalizations 
of the proposed DDW state. 

The present analysis is incomplete in one important 
respect. I have been unable to make firm statements re- 
garding the quasiparticle excitations and doping in the 
present model, although it is clear that power-law bond 
current correlations in the disordered phase probably 
leads to non-fermi liquid behavior of the electronic ex- 
citations. Since, we have considered models of the spin- 
singlet variety, there cannot be any spin-charge separa- 
tion by definition. The non-fermi liquid behavior must 
then be akin to that of the spinless Luttinger model in 



one dimension. One should note that a rigorous descrip- 
tion of electronic excitations are also unavailable for the 
RVB quantum dimer model. 

In the future, it may also be interesting to allow for 
sources and sinks, as they may represent an effective 
means of incorporating decay of collective excitations in- 
volving the order parameter into low Jxing nodal quasi- 
particles in the actual physical systeirEa. From this per- 
spective, the correct quantum model is actually a 8- 
vertex model in a transverse field, as discussed above, 
but without any complications of the projection opera- 
tors. 
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